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Extended Josephson junction driven by dc bias current is studied numerically. Two types of solutions, symmetrical and asymmet- 
Q^rical, are found. The current-voltage characteristic (IVC) is calculated. The symmetrical solutions form main histeretic IVC and 
asymmetrical one create an additional branch. Depending on the bias current value periodic, quasiperiodic and chaotic modes of 
the junction motion was observed. Dynamics of the junction affected by thermal fluctuations was analyzed. Stability of different 
states of the junction is discussed. 
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2. Mathematical model 

To model the EJJ the two-dimensional (2D) damped sine- 
Gordon model was used |l|] 
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Extended Josephson junction (EJJ) driven by dc current is, 

S from the theoretical standpoint, a nonlinear dynamical system. 
Mathematical model to describe the EJJ is the damped sine- 
Gordon equation Depending on the current value the EJJ 
can be in a statical or a dynamical state. In the wide range of 
the damped parameter the junction IVC is hysteretic. In this 
_ case at some values of bias current the EJJ can be in one of two 
different states. In the absence of external perturbations each of 
~~ * these states is stable. But thermal fluctuations can switch the 
junction from one state to another fl. 

—i Stability of the junction states under thermal fluctuations is 
^vQ special interest in last years. Authors of the work have in- 
vestigated the rate of thermally induced escape from the zero- 
voltage state in long Josephson junctions (LJJ). Return to zero- 
voltage state current statistics under thermal fluctuation was an- 
alyzed in [4]. Lifetime of superconducting states of long and 
short JJs was investigated in works JD, 0] . 

The aim of this work is detailed study of the dc current driven 
EJJ in the in-line geometry without external magnetic field. 
Symmetrical and asymmetrical solutions of corresponding sine- 
Gordon equation were found numerically. The hysteretic IVC 
of the junction formed by solutions found was calculated. The 
asymmetrical solutions form the additional branch of the IVC. 
Depending on the dc bias current value the EJJ can demonstrate 
different types of motion. Periodic, quasiperiodic and chaotic 
behavior of the junction were revealed. The effect of thermal 
fluctuations on the stability of symmetrical and asymmetrical 
states was investigated. 
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Here ip is the phase difference of order parameters of the su- 
perconductors forming the junctions, j c is the critical current 
density of the Josephson junction, h is the Planck constant, e 
is the value of the electron charge, c is the light speed, C and 
g are the capacity and the conductivity of the junction per unit 
aria, d is the thickness of the junction. The random force F 
conjugated with the phase difference simulates effect of ther- 
mal fluctuations on the junction. It is the Gaussian ^-correlated 
random variable 



(F(r,t)F(r',t')) 
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In this work for the sake of simplicity the model junction in 
the form of an infinite L-wide strip is considered. The strip 
is directed along the y coordinate. The edges of the strip are 
defined as x — +L/2. The boundary conditions in the case of a 
dc current driven junction in the in-line geometry are given by 



dtp I 
toT 



=§.*') 



4n 2 d 



(4) 



Email address: artemovOfti.dn.ua (A. N. Artemov) 



where <E>o is the magnetic flux quantum and / is the value of 
the dc bias current per unit length of the junction. In the y di- 
rection the periodic boundary conditions with the period L are 
supposed. 
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To obtain dimensionless expressions r = (hC/2ej c ) 1 ^ 2 was 
used as the time unit, Aj - ^<&Qcl{%n 2 dj c ) was used as the 
length unit, voltage is measured in 2e/hr and current was nor- 
malized by AjLj c . In order to solve the equation with the ran- 
dom force numerically it was rewritten in the form of two dif- 
ferential equations of the first order in time 



(fi, — u, 

u, + yu + sin ip = tp xx + ip yy + frfy 
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with the boundary conditions 

<p x (±L/2,y,t) = ±1/2, 
u(±L/2,y, t) = <p t (±L/2,y, t), 



<p(x, 0, t) = ip(x, L, t) (7) 
u(x, 0, f) = u(x, L, f). (8) 



Here u(r, t) is the dimensionless local voltage on the junction, 
y is the dimensionless damping parameter, fo is the Gaussian 
random variable with the unit dispersion and the dimensionless 
dispersion is given by the expression 
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where Hj = Qo/nAjd is the overheating field of the EJJ, 
r = T/T c , T c is the superconducting transition temperature of 
the metal forming the junction. To extract the temperature fac- 
tor in eq.® the model temperature dependence of the London 

r2 was supposed. The Ginzburg- 



length A L (T) = A L (Q)/ VT" 
Levanyuk parameter 

urn 

can be written in the form 
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It characterizes the dependence of the random force simulating 
thermal fluctuations on material parameters of a junction. To 
gain some insight of this parameter value we evaluated it for the 
lead junction. Taking into account that in this case T c = 7.2K 
and A L (0) = 37nm we find Gi » 2 ■ 1(T 9 . 

In the absence of fluctuations the model considered is identi- 
cal to the LJJ because its solutions don't depend on the coordi- 
nate y. But the dispersion <£3j due to its dimension corresponds 
to the 2D space and ID time distributed random force. To con- 
sider a long (L » Aj) and narrow (w < Aj) junction one has 
to average eq.© over the width. Corresponding dispersion de- 
pends on w. But in this case in the in-line geometry we have to 
use more complicated boundary conditions. 

To solve eqs.d7][8]l the explicit finite difference scheme of the 
first order in time and the second order in space was used. The 
time interval was accepted At = 0.01 and the space interval was 
Ax = 0.1. 

3. Main branches of IVC without fluctuations 

In this work we intend to discuss the influence of thermal 
fluctuations on a behavior of the EJJ. Intensity of the fluctu- 
ations depends on temperature. In our model this fact is taken 
into account as the temperature dependence of the random force 
dispersion (|9). In this situation the effective junction width 
L depends on temperature because it is measured in the unit 



Aj(T). The time independent width of the junction considered 
is Lo = 14.14/1/(0). At temperature T = 0.5T C which is ac- 
cepted in this work the effective width is L = 10. 

At first, we examine the EJJ behavior in the absence of fluc- 
tuations. In this case solutions of the EJJ don't depend on the y 
coordinate and in fact the junction is long. To get overview 
about dynamical states of the contact we have calculated its 
IVC. We consider the IVC as a diagram in the / - (u) plane, 
where (u) is the dc voltage, each point of which depicts the 
attractor of the dynamical system. The diagram allows us to 
systemize the attractors with regard to a way in which it may be 
reached by the system. 

To calculate IVC we start from the zero initial conditions at 
current value less then critical one, which is I c = 4. Then we 
increased it by small steps up to maximum value to obtain the 
ascending branch of the IVC. In a similar way starting from the 
attractor at highest current value and decreasing current step by 
step up to reaching a zero voltage state, the descending branch 
was obtained. Solutions at the previously calculated nearby 
point are taken as the initial conditions for the next point. After 
each changing of the current the system during t = 500 was re- 
laxed to a new attractor. Then the dc voltage (u) on the junction 
was calculated as the averaging of the value of the local voltage 
u(r, f) over Lx L space area and t = 1000. 
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Figure 1: IVC of the EJJ. Arrows show the the direction of the current varying 
on the IVC branches. Dashed line plots the separated branch of the IVC. Insert 
shows the portion of the IVC confined by the dashed rectangle. The segments 
of the branches, where quasiperiodic and chaotic types of motion are observed, 
are pointed. 

Resulted IVC is shown in figQjby closed circles. As it fol- 
lows from the way in which the IVC was obtained, all these 
attractors can be reached starting from zero initial conditions. 
The IVC demonstrates hysteresis behavior. A typical peculiar- 
ity of such IVC is that the return into zero voltage state current 
is less than critical one. For EJJ discussed it is / » 3.837. 

Solutions of eqs.(|5]|6]l positioned on the main IVC possess 
the space symmetry specified by the zero initial conditions and 
the boundary conditions (0|8]l. That is the phase difference and 
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the voltage distributions are symmetrical with respect to the 
center of the junction (x = 0) and the magnetic field distribution 
is antisymmetrical. Such solution are well known QjJ, lS|] . 

To classify the type of the dynamical system motion at the 
given bias current / the phase portraits u(r, t) vs cos cp(r, t) was 
plotted, were r is the space coordinate positioned in the center 
of the junction. In the most points of the IVC the motion of the 
system is periodic. The attractors of such type of the motion 
are closed lines. We will not discuss this type of motion in 
more details here. 

To analyze more complicated motion we calculated the time 
series v, = (u(r, f,)) at given current /. Each term of the series 
is the voltage averaged over L x L junction area calculated with 
the time interval t = 1. Each such series contains 5000 terms. 

Complicated multi-frequency motion of the system is ob- 
served in the current interval / e (4.62-5.26) on the descending 
branch of the IVC. The trajectory on the phase portrait visually 
fills up the whole region. The power spectrum of the voltage 
for such motion is shown in fig|2]and it demonstrates multitudi- 
nous spectrum lines. The evaluation of the correlation dimen- 
sion gives £>2 = 2 what is an argument for a quasiperiodical 
motion. Such a motion consists of periodic motions with in- 
commensurate frequencies. 
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Figure 2: Power spectrum of the quasiperiodic motion at the bias current value 
7 = 5.0. 

Rather complicated picture is observed in the current inter- 
val from / = 3.881 to the return current / = 3.837. Here nar- 
row current intervals where the system demonstrates a chaotic 
behavior are intermited by those with quasiperiodic motion. 
Thus, at currents / = 3.881 - 3.865,3.854,3.849,3.844,3.840 
the chaotic motion was found, whereas at / = 3.864 - 
3.855,3.852,3.847,3.842 the quasiperiodic one is observed. 
The power spectrum of the chaotic motion at current / = 3.87 
is shown in fig|3] It is characterized by the wide band of a con- 
tinuous spectrum what is typical for a chaos. The correlation 
dimension of the strange attractor is evaluated to be D2 — 1 .73. 

All attractors placed on the main IVC are joined by the fact 




Figure 3: Power spectrum of the chaotic motion at the bias current value / = 
3.87. 

that they can be reached starting from symmetrical zero ini- 
tial conditions. These conditions are similar to those which are 
take place in real experiments with EJJ driven by dc bias cur- 
rent. And this IVC contains all possible dynamical states of the 
junction which can be reached by means of changing bias cur- 
rent without additional external influences. But the dynamical 
system (|5]|6]l has another attractors which don't contain sym- 
metrical initial conditions in their domains of attraction. Some 
of them will be discussed in the next section. 

4. Separated branch of IVC 

In previous section the main branches of IVC were consid- 
ered. Here an additional separated set of attractors are dis- 
cussed. This branch is shown in fig[T]by dashed line. It is 
placed in the low bias current (/ € (1.7 - 4.16)) region of the 
IVC. Nonzero voltage states can be stable at current much less 
than critical on the branch. The system can be moved reversibly 
along the whole branch by means of bias current varying ex- 
cluding the highest and the lowest current points where the sys- 
tem switches to the main branch. 

Solutions of eqs.((5]|6]) on this branch are asymmetrical in 
contrast to that on the main IVC. Domains of attraction of 
these attractors don't contain symmetrical (zero) initial condi- 
tions. Asymmetrical initial conditions must be used to obtain 
such solutions. Thus, the initial phase difference distribution 
(p{x) = 0.2x 2 + 0.6x + 0.4 and zero voltage distribution fall 
into the domain of attraction of the attractors on the separated 
branch corresponding to the current interval / e (3.7 - 4.0). 

The solutions in the upper segment of the branch (/ > 4) 
correspond to states in which two fluxons enter in the junction 
from opposite edges in different moments of time and annihilate 
inside of the junction. On the lower segment one only fluxon 
is situated in the junction. It enters from left edge, moves to 
right and exits (annihilates) from right edge. Then the fluxon 
of opposite sign enters through right edge and moves to left. 
Space-time evolution of such type of solution is shown in fig@] 
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On the most part of the branch a periodic motion is observed. 
In this case the period of motion consists of two stages of en- 
tering and annihilation of fluxons. 

The system demonstrates also a chaotic behavior in the cur- 
rent interval / e (3.98 -4. 10). The power spectrum correspond- 
ing to the bias current value / = 4.02 is shown in fig|5] The cor- 
relation dimension of the attractor at such current is Di — 1 .86. 




Figure 4: Magnetic field space-time evolution at / = 3.85 on the separated 
branch of the IVC. 
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Figure 5: Power spectrum of the chaotic motion at the bias current value / = 
4.02 on the separated branch. 

Existence of the separated branch of attractors makes it pos- 
sible for the EJJ to demonstrate different types of motion at the 
same value of current. We focus your attention on the fact that 
in the current interval / e (3.98 - 4.10) the system demon- 
strates periodic type of motion on the main branch of the IVC 
and chaotic one on the separated that. 



5. Influence of thermal fluctuations on the JJ behavior 

Usually the influence of thermal fluctuations on the dynam- 
ical behavior of a system depends on the ratio between typical 
sizes of an attractor and fluctuations. If the size of fluctuations 
is small in comparison with the attractor one the dynamical mo- 
tion can be separated from noise and analyzed. In the opposite 
case fluctuations destroy the attractor and noise will predomi- 
nate. 

To investigate the influence of thermal fluctuations on the dy- 
namical properties of the EJJ we solved eqs.((5]|6} numerically 
taking into account the random force. The intensity of model 
fluctuations was defined by the parameter Gi = 10~ 7 and the 
reduced temperature t = 0.5. 

Results of the calculations can be summarized as the follow- 
ing. Characteristic sizes of the fluctuations under conditions 
given are noticeably less than attractor ones. So, the system 
moves in the vicinity of the attractor corresponding to given 
bias current. Dc voltage on the IVC branches remains practi- 
cally the same as without fluctuations. 

But presence of the separated set of attractors made the sys- 
tem behavior to be very sensitive to the fluctuations influence in 
the bias current interval / e (3.84 - 4.10) where two solutions 
of eqs.(|5]|6]), symmetrical and asymmetrical, can be realized at 
each value of current. It is found that under given conditions the 
lifetime of symmetrical solutions (on the main IVC) are much 
less than that on the separated branch. The thermal fluctua- 
tions are large enough to push out the system during finite time 
(t ~ 10 3 ) from a position on the main branch into the domain 
of attraction of the attractor on the separated one. At higher 
currents (/ > 4.1) the stability of the system states on the main 
IVC rises steeply. Switches of the states was not observed dur- 
ing t = 2 ■ 10 s . 

One of the possible variants of actual IVC is shown in figj6] 
by lines marked by open circles. It was calculated in the same 
way as before. The picture represents the reduced portion of 
the IVC on which difference between actual and fluctuationless 
IVC is visible. At more high values of the bias current the IVC 
practically doesn't depend on fluctuations. 

The process of the system switching from one metastable dy- 
namical state to another is probabilistic. It means that actual 
IVC will depend on the regime of calculations or measuring. 
During calculations under regime chosen (f = 500 + 1000 for 
each point) we don't observe revers system hops, from the sep- 
arated into the main branch. The more detailed examination of 
the states stability (calculation time t ~ 5 • 10 s ) doesn't estab- 
lish revers hops too. It means that the lifetime of the system 
on the separated branch is t s » 10 5 whereas that at main one 
is T m ~ 10 3 - 10 4 . Switches between metastable states on the 
ascending and the descending branches of the main IVC also 
was not observed. 

6. Conclusion 

Nonlinear systems at given value of driving parameter can 
have several attractors separated by separatrix in a phase space. 
To locate the system on a partial attractor it is need to specify 
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Figure 6: Dependence of I vs (u) for EJJ taking into account thermal fluctua- 
tions. Open circles mark ascending and descending branch of the actual IVC. 
Dashed lines plot the branches of fluctuationless main IVC. 

initial conditions within its domain of attraction which often is 
unknown. 

Well known sample of such a system is the Josephson junc- 
tion with hysteretic IVC discussed above. At some values of 
the bias current the system can occupy one of two attractors. 
A state of such a junction at given current depends on the way 
to obtain this state. All solutions on the main IVC possess the 
symmetry specified by the antisymmetric al boundary and the 
symmetrical initial conditions. 

In this work the additional set of attractors is found. It termed 
as the separated branch of the IVC. Corresponding solutions 
are asymmetrical. To find them asymmetrical initial conditions, 
similar to one-fluxon solutions, must be used. 

In the absence of external perturbation the solutions on both 
the main and the separated brunches are stable. Under fluc- 
tuations affect the situation is different. If at given bias cur- 
rent value the system allows more than one solution, they be- 
come unstable for switching between them. It is found that 
under conditions chosen (Gi = 10~ 7 , T - Q.5T C ) the lifetime 
of the system on the main branch at current / < 4.1 is quite 
small T m ~ 10 3 in comparison with that on the separated one 
t s » 5 • 10 5 . That is, in these conditions the system state on the 
separated brunch is preferable. At current / > 4.1 switching be- 
tween branches was not observed during time period t — 2 ■ 10 5 . 

Switching between the separated and the main branchs re- 
veals itself in the change of values of dc voltage and return cur- 
rent. Moreover, as a result of such a switch a type of motion can 
be changed. Thus, in the bias current interval / e (3.99 - 4.10) 
the EJJ switches from periodic attractor on the main IVC to 
chaotic that on the separated brunch. 

The IVC considered (figHJl is determined by the characteris- 
tic time (f = 500 + 1000) spent to calculate one point on the 
IVC. This time is rather large as a calculation one but it is much 
less then a characteristic time of a dc current experiment. The 
behavior of a real junction will depend on the fluctuations inten- 



sity (parameter Gi, temperature) and details of the experiment. 
It can vary from a steady dynamical state on the branch of the 
main IVC to a telegraph noise signal with two or even three 
levels of voltage. 
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